Quantized charge pumping by surface acoustic waves in ballistic quasi-lD channels 



Vyacheslavs Kashcheyevs, 1, Q Amnon Aharony, 1,2 and Ora Entin-Wohlman 1 ' 2,3 

1 School of Physics and Astronomy, Raymond and Beverly Sackler faculty of Exact Sciences, 
Tel Aviv University, Tel Aviv 69978, Israel 
2 Department of Physics, Ben Gurion University, Beer Sheva 84-105, Israel 
s Albert Einstein Minerva Center for Theoretical Physics 
at the Weizmann Institute of Science, Rehovot 76100, Israel 
(Dated: February 24, 2004) 

Adiabatic pumping of electrons induced by surface acoustic waves (SAWs) in a ballistic quasi- ID 
quantum channel is considered using an exactly solvable tight-binding model for non-interacting elec- 
trons. The single-electron degrees of freedom, responsible for acoustoelectric current quantization, 
are related to the transmission resonances. We study the influence of experimentally controllable 
parameters (SAW power, gate voltage, source-drain bias, amplitude and phase of a secondary SAW 
beam) on the plateau-like structure of the acoustoelectric current. The results are consistent with 
existing experimental observations. 
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I. INTRODUCTION 

Single electron transport through low-dimensional 
mesoscopic structures, driven by surface acoustic waves 
(SAWs), is a subject of active experimental El II 01, El E , 
IE 0, H I! and theoretical El El El El El [ll El [Hp 
research, with potential applications in metrology [lE 
and new computation technologies |2C| . In a typical ex- 
perimental setup, a quasi-one dimensional ballistic chan- 
nel is defined in a AlGaAs/GaAs heterostructure and a 
SAW is launched in the longitudinal direction at a fre- 
quency u/2-k of several GHz. Under appropriate con- 
ditions, the acoustoelectric dc current I exhibits a stair- 
case plateau-like structure as function of the gate voltage 
(which controls the depletion of the channel) and of the 
SAW power. At the plateaus, the current saturates at 
quantized values I = e(w/27r)m, corresponding to the 
transfer of an integer number m of electrons per each pe- 
riod of the SAW (here e is the electron charge) . The first 
plateau is the most flat and robust to changes in the con- 
trol parameters; the higher plateaus become less and less 
pronounced as the plateau number m is increased. In 
addition, the effect of factors such as source-drain bias 
H temperature El 01) gate geometry 0], a sec- 
ondary SAW beam H Q , and perpendicular magnetic 
field [3 on the staircase structure and the quality of the 
first plateau have been studied experimentally. 

In the experiment, the plateaus are observed below 
the conductance pinch-off, when electrons in the source 
and in the drain reservoirs are separated by a potential 
barrier. This observation forms the basis for the simple 
qualitative explanation of the quantized transport which 
has been proposed in the first experimental report EI 
and further refined in Refs. H El and El They argue 
that when the wavelength A of the SAW is comparable 
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with the size of the depleted region L (as it is in the 
experiments 0, H 0), a single potential well forms on 
top of the static barrier. This potential well then acts 
as a dynamic quantum dot, which can hold an integer 
number of electrons due to the Coulomb blockade effect. 
The captured electrons are transferred from one side of 
the barrier to the other , with p ossible quantization errors 
due to back-tunnelling El |ll| . In the above description, 
the formation of the quantum dot and the transport of 
the localized electrons are treated separately. Particular 
effects which have been studied theoretically within this 
picture are the non- adiabatic effects at the quantum dot's 
formation stage E3j an d the classica l dy namics of the 
already confined interacting electrons 

A different perspective on the problem has been sug- 
gested in Refs. Il9 and This approach relates the 
acoustoelectric transport to adiabatic quantum pump- 
ing of non-interacting electrons. The external poten- 
tial, generated by the SAWs and by the control gates, 
is viewed as a perturbation acting on a coherent quan- 
tum wire • The resulting "staircase" structure of the 
acoustoelectric current and its dependence on model pa- 
rameters within this approach have been studied in Refs. 
Eland El using the adiabatic approximation in con- 
junction with an exactly solvable one-dimensional (ID) 
tight-binding model. This theory yields a crossover from 
a non-quantized acoustoelectric transport to the quan- 
tized limit as the SAW power and/or the static barrier 
height are increased. Although this picture requires the 
Coulomb interaction in order to set the energy scale of 
the problem [111 Il2l Il5| . the main qualitative features 
of the experiment can be reproduced within a model of 
non-interacting spinless electrons 0, ^| . 

In this paper we extend the results of Ref. 0. The 
mechanism of quantized transport is elucidated by u sing 
a resonance approximation for adiabatic pumping |23j. 
Both current quantization and transmission resonances 
are determined by the quasi-bound states of the elec- 
trons captured by a moving potential well. New effects, 
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including the influence of a counter propagating SAW, 
static potential asymmetry and source-drain bias on the 
number and shape of the quantization steps, are consid- 
ered. We compare our qualitative conclusions with the 
published theoretical and experimental results. In par- 
ticular, tuning the amplitude and the phase of a weak 
secondary SAW is found to improve the quantization by 
accordance with an earlier experimental report Q. For 
this effect, we propose a new quantitative relation be- 
tween the phase and the amplitude of the optimal sec- 
ondary SAW which can be easily checked using existing 
experimental setups. 

The results are presented as follows. In Section [H] we 
describe the model and the algorithms for calculat- 
ing the adiabatic current. In Section IlIII we explain the 
formation of the integer plateaus and make quantitative 
analytic estimates by applying the resonance approxima- 
tion 

El 

to the model of Ref. E Building on these 
results we analyze in Sec. IIVI additional factors, not de- 
scribed previously, such as reflected SAWs, source-drain 
bias and gradual screening of the pumping potential. Fi- 
nally, a discussion of our results in the context of related 
work is presented in Sec. [V] together with several con- 
clusions. 



II. THE MODEL 
A. The tight-binding Hamiltonian 

The model consists of a nanostructure and two ID 
leads connecting its ends to the electronic reservoirs 
(Fig. . The leads are modelled by two chains of sites 
with vanishing on-site energies and nearest-neighbors 
hopping amplitudes — J. An electron moving in the lead 
has the energy E(k) — — 2Jcosfca, where k is the wave 
vector and a is the lattice constant. 

The nanostructure represents the SAW-affected part 
of the quasi- ID channel. It is described by N tight- 
binding sites n = 1, 2, ... N, that form a chain of length 
L = (N — l)a. The non-zero elements of its Hamiltonian 
matrix, H.o(t), are constant nearest- neighbors hopping 
amplitudes — Jd and diagonal time-dependent on-site en- 
ergies e n (t). The connection between the ideal leads and 
the time-dependent part of the constriction is introduced 
through a hopping amplitude — Ji {—J r ) between the left 
(right) lead and the site 1 (N) of the nanostructure. The 
resulting full Hamiltonian of the quantum wire will be 
denoted H. In the trivial case where J l — J r = J d = J 
and e n = 0, H describes an ideal ID tight-binding lat- 
tice. The time-dependence of the on-site energies e„(t) 
reflects the SAW-induced piezoelectric potential, which 
is responsible for the pumping effect. The simplest case 
is that of a single running wave as introduced in Ref. llH 

e n (t) = -V g + P cos (u>t - qx n ) . (1) 

Here, V g is the gate voltage (in energy units), q is the 
wave-vector of the SAW, and P is the amplitude of the 




FIG. 1: One-dimensional tight binding model for SAW- 
induced pumping. 



piezoelectric potential, induced by a SAW running from 
left to right (for q > 0) . The origin of the spatial coordi- 
nate x n = (n— l)a — L/2 is chosen to be at the middle 
of the channel. 



B. The acoustoelectric current 

The adiabatically-pumped current flowing between two 
reservoirs with equal electrochemical potentials IM = Ur , 
is usually calculated using the Brouwer formula [24. l2a . 
We use an equivalent formalism, developed in Ref. l22l 
which also includes the effects of a finite bias eVsD = 
Mi — Mr 7^ 0. The total instantaneous current, I a (t), of 
spinlcss electrons from the lead a = l,r into the nanos- 
tructure consists of two parts, I a (t) = /P um P + 7^ las . 
These two parts can be conveniently written down using 
the instantaneous scattering states |\& Q (t)) (normalized 
to a unit flux), the instantaneous transmission coefficient 
T(i) and the overall scattering phase 6(t) of the nanos- 
tructure, [23 

C mp = ^j dE (*"(*)! « !*«(*)> d(fl d + E fr) ' (2) 

(3) 

Here f a (E) = 1/[1 + e 0( E -^ [ s t h e Fermi distribution 
with j3 = 1/fcsT (T is the temperature). If the system 
is unbiased, then i^ las = and Eq. Q can be shown 
to reproduce [12, H^l the Brouwer formula HI • On 
the other extreme, if no pumping potential is applied, 
/P ump = and Eq. © leads to the Landauer formula 
pfij for the conductance, G = (e 2 /h)T. 

For most of the discussion we assume both the bias 
voltage Vsd and the temperature T to be zero. In this 
case only electrons at the Fermi energy /J-i — fi r = Ep 
participate in the scattering. Solving the scattering prob- 
lem for the potential and using Eq. @ yields the 
charge Q pumped over one period (the average dc com- 
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ponent of the current) |l6j . 

eJrsinka f 27r/ " N 



Q 



dtif nmv (t) 



sinfca r*'" ^ . . 

/ dt },£n\gn,l\ 



(4) 

[g- 1 ] n n , = [EI - H Q } ntn , + 8 n<n , e lka (* n>1 J, + 5 n , N J r ) , 

(5) 

where J; ;r = Jf. r / J and fc is the Fermi wavenumber, 
i?p = E{k). The instantaneous transmission is 

T(i) = 4| 5JV ,i| 2 J;J r sin 2 fca. (6) 

The integrand in Eq. Q is a meromorphic function of 
z = exp(iu>t), with 2N pairs of complex conjugate poles. 
Therefore, the integration of Eq. ma y be carried out 
exactly, once the positions of the poles are determined by 
solving numerically the corresponding algebraic equation 
of degree 2N. 



C. Resonance approximation 

The second term on the r.h.s. of Eq. (JSJ is the self- 
energy addition to the Green's function of the isolated 
channel, due to the coupling to the external leads. When 
the latter is sufficiently small, the total pumped charge 
can be divided into contributions from separate single- 
particle levels of H.Q. A systematic development of this 
approach leads to the resonance approximation for pump- 
ing, which is discussed in detail in Ref. [23l Here we 
summarize the resulting algorithm for calculating the 
pumped charge in the this approximation. 

1. Solve the instantaneous eigenvalue problem 



(m) 



E m ip 



(m) 



and obtain the 



approximate resonance energies E m (t). 

2. Calculate the time-dependent decay widths of each 
resonance into each lead, 



| r M,rM} = {j i |^ T ' 



sin ka . 



(7) 



3. For each to, find all such times t m j at which the 
resonance condition E m (t m> j) — Ep is satisfied. 

4. At each resonance time t — t m j, compute the par- 
tial charge transferred between the left lead and the 
mth quasibound state in the channel, 



AQ r 



(m) 



-,(m) 



-,(m) 



(8) 



5. Calculate the total charge pumped from left to 
right 1 : 



= - ^ AQ mj - sgn E m (t m ,j) ■ 



(9) 



or set 
step 3. 



if no resonances were found in 



The algorithm has a direct physical interpretation 23] . 
Whenever the energy E m of a (quasi-)bound state crosses 
the Fermi level Ep, an electron either occupies ("load- 
ing") or leaves ("unloading") this state. The correspond- 
ing unit pulse of current is distributed between the chan- 
nels proportionally to the ri^'s. Except for specifically 
designed Hamiltonians Ho(t), Q ICS — > Q in the limit of 
vanishing couplings r„ — > 0. 

The resonance approximation fails when either (i) the 
total width of a particular resonance is larger than the 
distance to the next energy level; or (ii) the partial de- 
cay widths change considerably while the system is 
at resonance j23|. As discussed in detail in the follow- 
ing section, these restrictions become significant for the 
non-quantized transport, but have little influence on the 
shape of the current quantization steps. In all the cases 
in which the resonance approximation is inadequate, we 
rely on the results of an exact calculation. 



III. FORMATION OF QUANTIZATION STEPS 
A. Application of the resonance approximation 

The results of a full calculation (as outlined in 
Sec. IIIB|) show that the pumped charge, Q, fol- 
lows a staircase- type dependence on the gate voltage, V g , 
and/or on the SAW amplitude, P, for a wide range of the 
model parameters. This 'quantization' can be related to 
the structure of the transmission resonances [2^, H3, . 
We first establish this relation quantitatively and then 
use it to analyze various aspects of the model. 

The calculation of the pumping curve can be visual- 
ized using a diagram like the one shown in Fig. |5| First, 
one plots the instantaneous eigenvalues E m for V g = as 
function of time uit (curves in the right panel of Fig.|2J). 
The small circles on the top of each curve show the time 
evolution of the corresponding partial charge AQ: the 
diameter of each circle is proportional to \AQ/e\ < 1; 
shading is determined by the sign — black (•) for E m < 
("loading") and white (o) for E m > ("unloading") 
. Once the eigenvalue diagram is constructed, the set 
of resonances for each particular value V g of the gate 



1 Due to charge conservation, it is sufficient to calculate the charge 
transfer from the left reservoir. Therefore, the channel index a 
is fixed to a = / in Eqs. JHJ and . 
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FIG. 2: (Color online) Construction of the pumping curve Q Tes (V g ) in the resonance approximation. Right panel: instantaneous 
energy levels of Ho at V g = as function of time. Left panel: pumped charge Q res (in units of e) as function of gate voltage. 
Horizontal dashed lines show the correspondence between sharp features in the pumping curve (left), and the change in the 
number of energy levels at resonance (right); see text for a detailed discussion. All energies are given in units of Jd\ the 
parameters of the potential are: P = 8Jd, A = 4L, TV = 10. 



voltage is determined graphically: a horizontal line with 
ordinate Ep + V g crosses the eigenvalue curves in the 
right panel at the points where the resonance equation 
V g + E m (Vg = 0) = Ep is satisfied (step 3 of the algo- 
rithm). The abscissas of the crossing points determine 
the resonance times t m j to be used in Eqs. JHJ and 0. 
(The dashed horizontal lines in Fig. [21 mark the extrema 
of the eigenvalue curves, and thus correspond to par- 
ticular values of V g at which the number of resonances 
changes.) Finally, the total pumped charge, Q res (V^), is 
calculated by summing up the contributions to Eq. 10: 
the magnitude and the sign of each term is given by the 
small circle at the respective crossing point in the right 
panel. The resulting pumping curve Q rcs (V g ) is plotted 
in the left panel of Fig. [21 

Several aspects of the model are illustrated by the con- 
struction in Fig. [21 One can see that the quantization of 
the pumped charge is caused by electronic (hole) states 
with the lowest (highest) energy. When resonances oc- 
cur, (namely, at {t m ,j}), these states are localized near 
one of the channel exits — either Ti /T r <C 1 or Ti /T r ^S> 1 
— and therefore transfer almost integer charges [Eq. ©]. 
The number of steps counts the number of localized states 
involved. 

The exact integration [Eq. (@J] takes into account the 
"external" parameters of the model, ka, Ji and J r , which 
are ignored in the resonance approximation. In the 
following, we will consider only symmetric couplings, 
Ji = J r = J. We have calculated the exact pumped 
charge, Q(V g ), for several values of the "external" param- 
eters, but with the same pumping potential as in Fig. [21 
Representative results are shown in Fig. [21 along with the 
approximate Q Tes (V g ) from Fig. [21 [thin (blue) line]. For 
sufficiently small J, the exactly calculated curves contain 
integer steps and sharp, non-quantized features at large 
values of V g (e.g., the spikes marked by small arrows in 



Fig.OJ. The first steps are robust and do not change their 
positions as J and ka are varied (except for a trivial shift 
of Ep). The top of the pumping curve and the spikes are 
more vulnerable: as J is increased, the upper steps and 
the sharp features shift and become rounded. Narrow 
spikes disappear for J = Jd and ka close to the center of 
the band [see curve (c) in the right panel of Fig. [3]. 

The resonance approximation reproduces all the de- 
tails of the exact calculation for J <C Jd, because the 
resonance widths in Eq. J7J vanish in the limit of J — > 0. 
The non-generic sharp features are determined by the 
surroundings of level anti-crossings (see Fig. [2J, where 
the corresponding level spacings are tiny. As we expect 
form the validity condition (i) in Sec. Ill CI the finite res- 
onance width effects arc most important in this region. 
Indeed, the discrepancies between the exact and the ap- 
proximate curves in Fig. are well correlated with the 
fact that the shifts and the widths of the resonance lev- 
els for a tight-binding model are proportional to J cos ka 
and Jsinfca, respectively [Eq. lf7|]. 

We have made a similar comparison between the exact 
integration and the resonance approximation for several 
sets of "internal" parameters, P/Jd, A/L, and N. The 
most important conclusion is that the stair-case struc- 
ture of the pumping curve can be reliably estimated us- 
ing the resonance approximation. Hence, we will utilize 
this useful technique in the following as a source for vari- 
ous analytic estimates that will be further checked versus 
exact calculations. 



B. SAW parameters and the number of 
quantization steps 

For the lowest part of the spectrum (which is relevant 
for the quantized transport), the on-site energies Q can 
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FIG. 3: (Color online) Pumped charge versus gate voltage for different external parameters: J = 0.16 Jd (left panel), J = Ud 
(right panel), calculated in the resonance approximation (a) and exactly for the bottom of the band (b) and at the band center 
(c). Parameters of the potential are the same as in Fig. [21 



be treated as a potential function of a continuous spatial 
coordinate x n . For A > L, only one minimum of this 
potential can be located inside the SAW-affected part of 
the channel. The position of this minimum xo(t) = (t — 
t(j)v moves with the sound velocity v = to/q and passes 
through the middle of the channel at time to = w _1 (7r + 
2ir x integer). Electronic states localized in this moving 
potential well can be approximated by simple harmonic 
oscillator wave- functions jlfij . The corresponding energy 
spectrum is E m = — P — 2Jd + A(to— 1/2), m = 1,2, . . ., 
with a constant spacing A/Jd~ qa\j2P/ Jd- The lowest 
energy wave-function is approximately a Gaussian, 

4> [ n ] {t) = (C 2 k/2)- 1/A exp{-[x n - Xo (t)] 2 /e}, (10) 

with £ = 2a\jJd/A. The localization length of the higher 
levels can be estimated as £ m = £ \pm. 

The harmonic approximation is valid as long as the 
wave-packet is driven adiabatically by a parabolic well 
and is not perturbed neither by the ends of the channel, 
nor by the "hills" of the cosine-shaped potential profile. 
This implies the validity condition 

£ m «min[£/2-|*o(t)|,A/2]. (11) 

In order to illustrate the above reasoning, we draw the 
attention of the reader to a set of constant and equidis- 
tant energy levels E m (t) in the right panel of Fig- El i n 
the vicinity of ut = it. The lowest energy level follows the 
harmonic approximation as long as the parabolic mini- 
mum is located inside the channel, that is for the fraction 
X/L = 1/4 of the full period. Higher energy levels remain 
constant for shorter times, since their respective localiza- 
tion lengths entering Eq. Ijlll) are longer. The harmonic 
structure of the energy levels translates into a sequence 
of equidistant steps in the pumping curve, Q(V g ), with 
the same energy spacing A, as shown in the left panel of 
Fig. At each value of gate voltage, Vg = E m — Ep,& 
new pair of resonances and another step in the pumping 



curve emerge. The plateaus are rather flat because the 
resonant states at the loading (unloading) moments are 
well localized at the entrance (exit) of the channel. 

The number of quantization steps, A s t ep s, is limited 
by two competing mechanisms. The first limit is set by 
the number, Ni, of localized states that can be trans- 
ferred below the Fermi energy. If x n can be considered 
as continuous, the localization condition is roughly the 
same as the validity condition (|ll|l for the harmonic ap- 
proximation. For L < A it follows from £jvi = L/2 
that Ni = L 2 A/(l6J d a 2 ) = (71V2/8) N {L / \) yfP/T d . 
On the other hand, for large enough P the discrete- 
ness of the tight-binding grid cannot be neglected. For 
a rough estimate, we assume that the continuous ap- 
proximation breaks down if it yields an average distance 
£ m /m between the successive zeros of the mth wave- 
function, which is smaller than the inter-site spacing 
a. This happens for m > N 2 , where N 2 = N 2 /(4iVi). 
Putting the two limits together we estimate the num- 
ber of quantization steps, -/V s t C p S , as the integer closest 
to min(iVi, 7V2). By adjusting the parameters one can 
obtain at best a sequence of N/2 steps. The o ptim al pa- 
rameters L = 0.3A, N — 6, P = 8J d of Ref. indeed 
yield N stcps w JVi = 2.83 w N 2 « N/2. The decrease 
in the number of steps with increasing L/X reported in 
Ref. corresponds to the tight-binding limited regime 

iVsteps « N 2 oc X/L. 

Despite a certain inherent uncertainty of our estimates, 
they prove useful for understanding the effect of changing 
the amplitude and the wavelength of the SAW (Fig. 
In Fig. the number of steps is close to optimal, 
N/2 = 5, and is limited by the localization criterion 
A^steps ~ Ni = 3.9. Increasing P by a factor of 4 (Fig. 0Ja) 
reduces the number of steps due to discrete lattice effects: 
A s teps ~ ^2 = 3.2. One can clearly see that for higher 
energy levels (close to the band center) the tight-binding 
coupling Jd is no longer relevant: E m {t) with m > N 2 fol- 
low a sequence of cosine curves P cos(ujt + 8(f) with equal 
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FIG. 4: Instantaneous eigenenergies E m (in units of Jd) for 
V g = 0, N = 10 and (a) P = 8J d , A = 4L; (b) P = 32J d , 
A = 4L; (c) P = 2J d , A = 4L; and (d) P = 2J d , A = 2L. 
Insets: the pumped charge Q res as function of the scaled gate 
voltage V g +EF; the distance between the ticks on the ordinate 
axis is equal to a unit charge. 



phase differences qa. These curves correspond to the indi- 
vidual on-site energies e n (t). In this regime the hopping 
amplitude Jd leads only to tiny anti-crossings between 
the energy levels, which in turn give rise to the sharp 
peak-like structure in the pumping curve. The effects of 
the tight-binding approximation become less pronounced 
as P is reduced below the optimum (Fig. ^J:). In this 
case -/Vgteps ~ Ni = 1.96 and the peaks in the pumping 
curve are suppressed. The missing steps can be brought 
back by shortening the wave-length, as shown in Fig.0Ji. 
The estimated number of steps is now the same as in the 
original case (a). However, the non-parabolic shape of 
the potential minimum is more pronounced. Note that 
in case (d) the flat region for Eo(t) extends over half of 
the period, since A = 2L. 



IV. PERTURBATIONS OF THE PUMPING 
POTENTIAL 

The pumping potential in Eq. Q is of a rather high 
symmetry. Small perturbations - such as a static impu- 
rity or a reflected SAW — can change the shape and the 
position of the current quantization steps. In order to 
explore these effects, we add to e n (t) a smooth function 
of x„ and t. 



e n (t) = -V g + Pcos(u>t - qx n ) + U(x n ,t) ■ 



(12) 



Similarly to the situation discussed above, the structure 
of the relevant energy levels can be analyzed using the 



harmonic approximation, provided that U(x n ,t) changes 
slowly and the travelling wave-packet is well localized: 
^ m dU(x,t)/dx <C A. The first-order approximation for 
the instantaneous energy, 

Em(t) = -P- 2J d + A(m - 1/2) + U(x (t), t) , (13) 



is valid for |xo(t)| 



< L- 



£ m (we consider the case A > L) 



Note that xo (t) is the position of the potential well min- 
imum, and Xo(to) — (the middle of the channel). Now 
even within the harmonic approximation E m (t) is ex- 
plicitly time-dependent and this time dependence maps 
onto the shape of the current quantization steps. To 
make a quantitative statement we note that the instanta- 
neous wave-function remains unperturbed in first-order; 
\tl>if(t)\ becomes greater than |'0i(OI at t = to- At this 
point, the partial decay widths are equal, r[ m ' = Tr > 
and the resonance approximation yields a half-integer 
pumped charge. Therefore, the transition between the 
consecutive plateaus takes place at the gate voltages 
y( m > = E m (to) — Ep. In particular, half of the first 
step in the pumping curve is reached at the gate voltage 



V 1/2 = V^ =V Q + U(0,h) , 



(14) 



such that Q(V 1/2 ) = e/2. Here V = -E F -P-2J d +A/2 
is the threshold voltage for the first step in the absence 
of perturbations. 

The resonance moment associated with the left-right 
transition at V g = Vi/ 2 is well defined, since the energy 
levels E m (t) are in general no longer constant in the vicin- 
ity of t = to. Therefore, the slope of the first quantization 
step can be estimated from the resonance approximation. 
The value of the total pumped charge at V g = Vys + 5V 
is dominated by the unloading resonance at t = to + St, 
where SV — Ei(to)5t. The other resonances contribute 
charges exponentially close to an integer; for simplicity, 
let us consider only one loading through the left lead 
(which gives AQi « e = const.) before unloading at to- 
The contribution of the latter, AQ2(t), can be calculated 
using the Gaussian wave- function (|l(Jfl in Eqs. J7J| and 
iJSJ. The resulting total pumped charge AQi + AQ 2 is 



Q 



\^\t)\< 



1 + tanh ■ 



LvSt 



We define the steepness of the first step, S, as 



S : 



dQ 

dVa 



v a =v, 



eLv 



dU 
dx 



dU 
~dt 



(15) 



(16) 



t=t 

x=0 



The pre-factor in Eq. (|16|l is the least accurate, since the 
applicability of Eq. (|10fl at the ends of the channel is 
marginal. Taking the absolute value in Eq. I|16|l makes 
the result valid for both signs of E\ at t = to. Our deriva- 
tion is not justified for perturbations that yield small val- 
ues of the denominator in Eq. (|16fl . Then, the steepness 
remains bounded due to the finite resonance width. 
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FIG. 5: Interference of two counter-propagating SAWs with 
an amplitude ratio P~ / P = 0.2. The phase difference tp is 
changed from to 7ir/6, in steps of 7r/6, as indicated. Solid 
curves show the pumped charge (in units of e) versus the 
gate voltage (in units of Jd), the dotted lines mark the best 
quantization conditions achieved at tp — and tp = n. For 
tp € [n, 2n] the pumping curves repeat the same sequence 
in reverse order (not shown). The parameters used are: 
P — 8Jd, A = 4L, TV = 10. Curves are computed using 
the resonance approximation. 
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FIG. 6: (Color online) The steepness of the first step S = 



dQ/dVg at 



e(l/2), in units of eJJ 1 for P~/P 



0.05, 0.1, 0.2 as function of the phase difference tp. The curve 
for P~ /P = 0.2 corresponds to the data presented in Fig. 
Inset: threshold voltage Vi/2(cp) — Vq versus tp, in units of 
Jd- Curves are computed exactly from Eq. @. Thick dashed 
(blue) lines show analytic estimates, given by Eqs. (I18H19H for 
the smallest amplitude ratio P~ / P = 0.05; the pre-factor N\ 
in Eq. 119H has been treated as a free fitting parameter. The 
parameters used are: P = 8Jd, A = 4L, TV = 10, ka = 7r/5, 
■7 = 1. 



The quantization accuracy can be estimated along sim- 
ilar lines. However, the results are less transparent since 
the energy levels involved are beyond the simple har- 
monic approximation. 

A. Sensitivity to the second SAW 

For a particular example of a perturbation which mim- 
ics the experimental situation, consider the following po- 
tential 

U(x n , t) = P~ cos (ujt + qx„ + ip) + bx„/L . (17) 

Here P~ and tp are the amplitude and the phase of a 
second SAW, propagating in the negative direction. It 
can be generated either due to reflections of the main 
beam |.2] or by a second transducer Q| . We also include a 
simple static perturbation [proportional to b in Eq. (|17fl ] 
which breaks the left-right symmetry of the channel in 
the absence of the SAW. The estimates in Eqs. (|14fl and 
(fT^f) become 

V 1/2 = V -P-cos<p, (18) 



\b + 2qLP~ sin^l ' 

[We have used the relation TVi = L 2 /( 4 £ 2 ) in GSM 

First we consider the case of a reflected wave only 
(6 = 0). A series of pumping curves for different val- 
ues of the phase difference is presented in Fig. [5] As can 
be seen from Eq. I)18|l. the threshold voltage changes pe- 
riodically in tp, reaching extremal values at tp = and tt. 
Between these special values of tp, the staircase structure 



is more smooth, and the steps are more symmetric: the 
convex and the concave parts of a step become almost 
congruent. The pumping curves are identical for ±tp due 
to the symmetry of the potential (|12l) with 6 = 0. 

For a quantitative characterization of the second SAW 
effect we have determined numerically the positions and 
the slope of the pumping curves at Q — e/2 without 
any approximations in Eq. fllf. The results are shown 
in Fig. |S| Tuning the phase difference tp for a fixed am- 
plitude ratio P^/P to the values at which the r.h.s. of 
Eq. I|19|) diverges enhances the steepness of the first step 
by orders of magnitude. The sharpest steps are achieved 
at the extrema of the threshold voltage V\/2, as shown in 
the inset in Fig. [S] and qualitatively in Fig. [SJ 

The above example shows that a symmetric pump- 
ing potential is favorable for quantization: the steep- 
est plateaus are achieved without a secondary SAW or 
with P~ 7^ and tp = 0, tt, when the total SAW po- 
tential V(x n ,t) = e n (t) is invariant under x n —x n , 
t — ► —t + const. 

Further reduction of symmetry is achieved by choosing 
6 7^ and P~ ^ 0. Here two regimes are possible. For 
small 6, the situation is similar to the previous case: the 
steepness is greatly enhanced at two values of tp between 
and 27T, when the denominator in Eq. I|19|) vanishes. 
In contrast, for b > 2qLP~ it is the static asymmetry 
of the channel that determines the slope of the steps, 
which now has only one wide maximum as function of 
tp. This behavior is illustrated in Fig.0 which shows the 
slope of the first step as a function of the second SAW 
amplitude and phase. The initial steepness at P~ = 0, 
6 = Ud is S = 12.2(e/Jd) for the selected model param- 
eters. That is more than an order of magnitude less than 
in the unperturbed (P~ = 6 = 0) case. Increasing the 
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FIG. 7: (Color online) The steepness of the first step for 
P = 2J d , A = 2L, N = 10, b = U d , as function of P~/P and 
tp. The dashed (blue) line marks the combinations of ampli- 
tude and phase at which the first order estimate [Eq. 1)19^ ] 
diverges. Note the logarithmic grey-coding scale. The steep- 
ness without the perturbation (P~ =6 = 0) is S = 613 ej Jd- 



amplitude of the second SAW improves the steepness for 
7T < tp < 2tt with a single wide maximum at tp ~ 3ir/2, in 
agreement with Eq. (|19(l . At P~ sa bj (2qL) the steepness 
peaks sharply, almost reaching the unperturbed value. 
Further increase of P~ reduces the steepness gradually, 
which now has two maxima in tp, approaching tp = tt 
and tp = 2tt for large P~, similarly to Fig. H3 This ex- 
ample shows that a weak counter-propagating SAW with 
properly chosen amplitude and phase is able to compen- 
sate for the static asymmetry of the pumping potential 
and significantly improves the sharpness of the current 
quantization steps. 

Available experimental data are consistent with our 
conclusions. Periodic oscillations of V1/2 have been ob- 
served in experiments with one active SAW transducer 
when the frequency of the SAW was varied. The period 
of these oscillations was found to match a full 2tt phase 
shift between the main SAW and a weak wave reflected 
from the other (inactive) transducer Later experi- 
ments, with two active transducers on both sides of the 
constriction, have confirmed this scenario, and a sequence 
of pumping curves similar to our Fig.[S]has been reported 
4] . Tuning of the second SAW amplitude and phase has 
enabled the authors of Ref. |4| to improve the flatness of 
the first quantization plateau. 

The key argument leading to Eq. I|19|) concerns the gate 
voltage dependence of the potential profile at the point 
where capture/release of an electron happens with equal 
probabilities from either side of the barrier. Therefore, 
the phase and amplitude dependence of the steepness, 
S(<p,P~), is expected to be insensitive to the particular 
choice of the pumping potential, as long as it leads to 
a clear sequence of current quantization steps. We sug- 
gest the following generic scenario of the plateau quality 
improvement, that can be checked by detailed measure- 



ments using existing experimental setups. One should 
measure the traces of the first step steepness S(tp) as 
function of the reflected wave phase tp for a set of gradu- 
ally increasing secondary beam amplitudes P~ . At small 
amplitudes, P~ < P~, the steepness is expected to have 
one broad maximum at some tp = tp Q . As P~ is in- 
creased, the value at the maximum, S(tpo), increases 
and at P~ = P~ , the maximum splits into two, S(tpi) 
and S(ifi2), with tp 12 — tp ± arccos(P t r/P _ ), as shown 
in Fig. [7| by the dashed line [for our model calculation 
tp a = 3tt/2 and P~ m b/(2qL)}. 

B. Source-drain bias and variations of screening 

Experimentally, acoustoelectric current can be studied 
along the full crossover, from the depleted to the trans- 
missive state of the quantum wire, by changing the volt- 
age on the depleting gate. Our discussion so far has been 
concentrated on the quantized single-electron transport, 
which is observed in the depleted regime. As the first 
conduction channel opens, the shape of the pumping po- 
tential in real space as well as screening effects become 
increasingly important pl| and the usefulness of our sim- 
plified ID spinless electron model is very limited. Keep- 
ing these limitations in mind, we will choose model pa- 
rameters that most closely correspond to a point contact 
near the depletion threshold, and illustrate the break- 
down of quantized transport. 

For P > Jd : = J, the tight-binding band is significantly 
deformed (see Fig. QJ, therefore we choose a relatively 
small SAW amplitude P = 0.5 Jd, but a large number of 
sites N = 24 to maintain N steps > 1. The Fermi wave 
number ka = tt/12 is taken close to the band bottom. 

Consider first the situation before the SAW is applied 
(P = 0). The zero-bias dc conductance of the chan- 
nel is determined by the transmission coefficient T (Lan- 
dauer formula, see Sec. Ill B(l . There is a potential barrier 
between the left and the right reservoirs for — V g > 0, 
therefore the value of V g + Ep = —2JdCoska w — 2Jd 
is expected to be the borderline between transmissive 
and blocked states of our channel. This corresponds to 
the depletion threshold of a true point contact. We plot 
the transmission coefficient in the absence of SAW ver- 
sus gate voltage, T(V g ), in Fig.|S^ with a thin (blue) line. 
For V g < 0, the transmission is exponentially blocked by 
a rectangular barrier of height « — V g and length L, while 
above the depletion threshold a Fabry-Perot-like pattern 
of high transmission is observed due to multiple reflec- 
tions at the sharp ends of the constriction. 

At a non-zero SAW amplitude, the transmission coef- 
ficient T(t) becomes time-dependent and the adiabatic 
formula [Eq. ©] should be used to relate it to the con- 
ductance. In the linear response regime, the second 
term in the curly brackets in Eq. © is proportional to 
d 2 f/dE 2 and can be neglect ed J22I . This results in a gen- 
eralized Landauer formula [22j, G = (e 2 /h)% iv , where 

T av = (oj/2tt) J Q dtT(t) is the time average of the in- 
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FIG. 8: (Color online) Crossover between the depleted and open regimes of the conduction channel for N = 24, ka — tt/12, 
J = Jd- (a) Time-averaged transmission coefficient %, v in the absence [thin (blue) line] and presence (thick black line) of a 
SAW with P = 0.5Jd, A = 2L. (b) The pumped charge Q [thick (red) line] and the total charge Q tot (thin black lines) for 
bias voltages eVsn/(tik>) changing from —3 to 3 in steps of one. The dashed-dotted line shows Q rcs at no bias [resonance 
approximation, Eq. JSJ]. (c,d) The same as in (a,b), but with an exponentially screened pumping potential [Eq. 121H 1. A 
sequence of snap-shots in panels (a,b) and (c,d) shows the corresponding pumping potential in real space for V g in the middle 
of the first plateau; time increases from top to bottom, uit £ [u)to — n /2uj;ujto + n/2]. 



stantaneous transmission coefficient T{t). This quantity 
is plotted in Fig. [SJl with a thick black line. One can 
see that switching on the SAW smears the sharp step in 
the conductance over the range of ±P around the de- 
pletion threshold. Qualitatively similar smoothing of the 
conductance quantization steps due to SAW has been 
observed experimentally 0, @] . 

Figure |Ht also shows some additional structure below 
the depletion threshold. This structure is correlated with 
the pumping curve shown by a thick (red) line in Fig. [SJd. 
Comparing Q(V g ) and T av (V g ) we see that each step in 
the acoustoelectric current is associated with a peak in 
the time averaged transmission as indicated by arrows in 
Fig. (the first two peaks are too small to be seen on a 
linear scale). It is easy to explain the origin of these peaks 
using the resonance approximation diagram (see Fig.[2Jl. 
At gate voltages between the quantization plateaus the 
system remains at resonant transmission for a consider- 
able fraction of the period, therefore T av becomes greatly 
enhanced. 

In the presence of both SAW and source-drain bias, 
the total charge transfer per period, Q tot = f^^dt Ii(t), 
becomes 

Q tot = Q + (e 2 VWM r av (20) 



(in the linear response regime). The result is a sum of 
the two terms: pure pumping contribution [thick (red) 
curve in Fig. |SJd]; and the average transmission (thick 
curve in Fig. multiplied by a constant. Equation 
(PU)l suggests that hu; is a natural unit for the source- 
drain energy mismatch eVsD- In the quantized pumping 
region, the contribution of the bias, J blas , becomes com- 
parable to that of pumping, _/"P um P ; if the bias voltage 
source transports several electrons per cycle. When T av 
is of order one, this regime is attained for eV equal to 
several tkj. Thus for a qualitative picture of the pump- 
ing curve in the presence of bias, we have plotted Q tot 
for the bias voltage eVso/tuj ranging from —3 to 3 by 
thin black lines in Fig.^D. The main observation is that 
the higher is the step number the more sensitive it is to 
the bias (as one can already appreciate form the aver- 
age transmission curve). Similar behavior is reported in 
experimental studies Q, ■ 

The main results of the above discussion remain un- 
changed if a phenomenological screening 0, |^| is intro- 
duced: 

e»(t) = [-V g + P cos (iot - qx n )} exp (-xl/L 2 s ) . (21) 

We have repeated the previous calculation using the same 
values of parameters but modified the pumping potential 
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(J2B with L s = i/4 = A/8. The results are shown in 
Figs. |U and (SJi. The main qualitative difference is the 
disappearance of the interference pattern in the transmis- 
sion curve both with and without the SAW. The number 
of steps is reduced to one (shown separately in Fig. |SJI'), 
since the effective amplitude of the SAW is decreased by 
the screening factor in Eq. i|21|) . Calculations with larger 
values of P produce more steps, along the same lines as 
discussed in Sec. IIII Bl for the unscreened potential QJ. 
We have also checked that the behavior of the first step 
steepness S(P~,ip) as function of reflected SAW ampli- 
tude and phase, follows the general scenario suggested in 
Sec. lTVAl 

We note that in the above example (Fig. |8J) the res- 
onance approximation still holds below the depletion 
threshold, when a moving quantum well is isolated from 
the Fermi sea in the leads. 



V. DISCUSSION AND CONCLUSIONS 

Quantized electronic transport, driven by SAW's, has 
been considered in several recent theoretical studies 0, 
ITU H3. IT3I IT3. ITU Il7| . Here we discuss our approach in 
relation to those works. 

Several of models 0, HTl IT2I Il5l | make a distinction 
between electrons already localized in a moving potential 
well (dynamic quantum dot) and those belonging to the 
Fermi sea. The current is then calculated by consider- 
ing the loss of electrons from the dynamic quantum dot 
at the stage of its formation 0, and/or its subse- 
quent motion 0, 0, 0, This approach presup- 
poses the formation of the dynamic quantum dot, but 
does not require it to be at thermodynamic equilibrium 
with the reservoirs at all times. Moreover, all the quan- 
tization error mechanisms within these models (gradual 
back-tunnelin g [ToL ITTL IT^| , non-adiabaticity at the for- 
mation stage Il2l and non-equilibrium dynamics during 
the transfer [jjl) consider electrons with energies that 
can significantly exceed the Fermi energy in the remote 
reservoirs. 

Our adiabatic quantum calculation [T^| differs from 
these studies in two significant aspects: (i) The forma- 
tion of a dynamic quantum dot is not a necessary con- 
dition for the calculation of the acoustoelectric current. 
We do identify, however, the localized electronic states 
(whenever such states are present) via the resonance ap- 
proximation and confirm that they are responsible for 
the quantized transport, (ii) In the adiabatic pumping 
approximation 22], the time-dependent potential never 
excites the carrier by a finite amount of energy away from 
the Fermi level |29| . Therefore, we never observe quanti- 
zation steps when the moving potential well rises above 
the Fermi level upon passing through the middle of the 
channel. 

We find the numerical calculations by Maksym [I3 to 
be the closest to our study. He considers a ID single- 
particle model with a pumping potential similar to our 



Eq. (|21() . The current at the quantization plateaus is 
found to be carried by the lowest energy states of the 
local potential minimum, in accordance with our results. 

The quantization accuracy in our approach is deter- 
mined by two factors which we expect to become experi- 
mentally relevant for sufficiently low tunnelling barriers. 
The first one is the possibility of both reservoirs to par- 
ticipate in the capture/release of an electron. This error 
mechanism is covered by the resonance approximation 
and leads to simple estimates like our formulas for the 
first step steepness discussed in Sec. IIV Al The second 
factor concerns mixing of the localized states with con- 
tinuum in the leads, which can give significant width to 
the quasi-stationary states in the moving quantum dot. 
Compared to the predictions of the loading/unloading 
scenario fSec. Ill C|) . this effect further degrades the flat- 
ness of the quantization (see, e.g., Figs. Eb, El an d EH) 
and eventually leads to the breakdown of the quantized 
transport as the channel opens fSec. HVB|) . 

We have not considered explicit Coulomb interactions 
between electrons in the depleted part of the channel, 
which set the energy scale of the problem. One can 
make a naive estimate of the level spacing A, which 
in the continuous limit is [ltl A = Kq^fPjmF . Using 
typical experimental values [j] for the SAW amplitude 
P = 20meV, wavelength A = 1/im and GaAs bulk effec- 
tive mass m* — 0.067mo, one gets A = lmeV, which is 
an order of magnitude less that the distance between the 
quantization steps observed in experiments This 
discrepancy can be qualitatively understood on a mean- 
field level: if an electron is captured by the moving po- 
tential minimum, its unscreened electric field makes the 
potential well seen by the other electrons much shallower, 
and thus increases the spacing A between resonances by 
the amount of the charging energy [23, 0OJ . Such a pic- 
ture is also supported by the numerical calculation of a 
two-electron problem by Gumbs and co-workers It 
is plausible that the effective values for the parameters 
of our model can be estimated from a self-consistent re- 
alistic calculation. 

In conclusion, we have considered a simple model for 
SAW-driven adiabatic pumping of electrons through a 
quasi-lD quantum wire. A stair-case structure of the 
acoustoelectric current has been mapped onto the instan- 
taneous energy spectrum of the pumping potential. Nu- 
merical calculations and analytic estimates confirm the 
experimentally observed behavior of the acoustoelectric 
current as function of the SAW amplitude, wavelength, 
source-drain bias, and the parameters of a weak counter- 
propagating beam. Quantitative measurements of the 
plateau quality as a function of the second SAW ampli- 
tude and phase are proposed to probe the relevance of our 
model. The presented single-electron picture captures all 
the main features of the quantized transport. 
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